Effect of bilayer coupling on tunneling conductance of double-layer high Tc cuprates 
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Physical effects of bilayer coupling on the tunneling spectroscopy of high Tc cuprates are investi- 
gated. The bilayer coupling separates the bonding and antibonding bands and leads to a splitting of 
the coherence peaks in the tunneling differential conductance. However, the coherence peak of the 
bonding band is strongly suppressed and broadened by the particle-hole asymmetry in the density 
of states and finite quasiparticle life-time, and is difficult to resolve by experiments. This gives a 
qualitative account why the bilayer splitting of the coherence peaks was not clearly observed in 
tunneling measurements of double-layer high-Tc oxides. 
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The interlayer coupling of electrons in high Tc cuprates 
was predicted to depend strongly on the in-plane momen- 
tum and vanish along the zone diagonals of the 2D Bril- 
louin zonei'^ This would lead to an anisotropic split- 
ting of the energy bands in bilayer compounds. This 
bilayer splitting was first observed by Feng et al in the 
angle-resolved photonemission spectroscopy (ARPES) 
of heavily overdoped (OD) Bi2Sr2CaCu208+:. (Bi2212) 
compounds.* The maximum splitting occurs near the 
anti-nodal points (±7r, 0) and (0, ±7r) and varies from 20 
meV in the superconducting state to 88 meV in the nor- 
mal state. Chuang et al also observed this splitting in OD 
Bi2212 samples in the normal state, but with a larger 
splitting energy 110 meV»^ Furthermore, by analyzing 
the energy dependence of ARPES spectra, Kordyuk et al 
concluded that the peak-dip-hump lineshape observed in 
ARPES are stemmed from the bilayer splitting.^ 

Tunneling spectroscopy is an important tool for ex- 
ploring low energy properties of high Tc superconductors 
(HTSC), as the tunneling conductance is proportional to 
the density of states (DOS) of electrons. The tunnel- 
ing measurements reveal important features of d-wave 
superconductors, such as the superconducting coherence 
peak at the gap edges and the V-shape low-energy spec- 
trum associated with the hnear D0S7i^i^i^°i^^i^^i^^ Be- 
sides these d-wave features, an asymmetric tunneling 
conductance background with a negative slope has also 
been observed iSifi In contrast to the ARPES, the bilayer 
coupling effect in tunneling experiments has not been re- 
ported. 

In this paper we investigate the effect of the bilayer 
coupling on tunneling measurements. The bilayer cou- 
pling splits the energy bands into the bonding and an- 
tibonding ones and leads to a separation of the super- 
conducting coherence peaks. This property can be used 
to probe the bilayer effect from tunneling measurements. 
Furthermore, it is shown that the DOS contributed from 
the bonding and antibonding bands behave differently in 
presence of the particle-hole asymmetry. With a negative 
slope in the normal density of states, the particle- hole 
asymmetry tends to reduce the DOS of the bonding band, 
but enhance that of the antibonding band. This will 
enhance one coherence peak but reduce the other one, 



and eventually cause the disappearance of two-coherence- 
peak structure. 

Our model involves the anisotropic c-axis coupling be- 
tween two Cu02 planes and is defined by the following 
Hamiltonian 



H = '^e(k)cl^^ci^i,^ + ^t_LCk){c\^j^^C2,ka + h.c.) 

Xl^W(4kT4-ki +Ci,_kiCi,kT), (1) 



i,ko 



where cj^^{l — 1,2) creates electons in the Z'th Cu02 
plane with momentum k and spin a. The kinetic en- 
ergy e(k) includes the chemical potential and thus the 
Fermi energy ep — 0. The superconducting energy 
gap is assumed to have d^2_y2 symmetry and A(k) — 
Ao(cosfca; — cosky)/2. In a tetragonal high Tc cuprate, 
the c-axis electron hopping integral is anisotropic >i*2i^ 

tj_(k) = —-^{coskx — cosky)'^. 

This anisotropy results from the hybridization between 
the bonding O 2p and unoccupied Cu 4s orbitals. If 
A(k) = in one of the double planes, Eq. ^ is the 
model that was widely used for studying the proximity 
effect in YBCO materialsi^ 

Defining the operators di^ko- = (ci^ko- + C2.k(T)/V2 and 
rf2,kcr ~ (ci,k(T ~ C2,kcr)/\/2, wc cau dccouplc the above 
Hamiltonian into two independent parts. 



H = H1+H2, 



(2) 



kcr 



-E^(k)(4kT4-ki + /»-c.), (3) 

k 

with ei^2(k) = e(k)±f_L(k). Hi 2 are the BCS Hamiltoni- 
ans for the bonding and antibonding bands, respectively. 
From the decoupled Hamiltonian, we can readily obtain 
the energy spectra of the Bogoliubov quasiparticles, 

H^ ^^Ei[\<i) (7-,kT7*,kT +7-,-ki7^.-ki) , 
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where Ei^2{'k-) = ^/A^(Xjll^(^Y±TJk))^ and the Bo- 
gohubov quasiparticle operators are defined by 



i,-ki' 



ki' 



and 



2 ^ i?,(k) 



(4) 
(5) 

(6) 
(7) 



The energy spectra of the two quasiparticle bands 
demonstrate the bilayer spUtting as observed in 
ARPESA^ Because of the anisotropy of t± (k) , the largest 
splitting occurs at the Fermi surface in the vicinity of 
(±7r,0),(0,±7r). 

The density of states of electrons p{llj) is defined by the 
imaginary part of the retarded Green function Gf\. (k, us) 

of fermions {d- ^cr' '^i.ko-}, 



p{uj) 



^5]lmG5(k,c^), 



ttN 



(8) 



where N is the total number of k vectors in the first 
Brillouin zone. For the i'th band, the density of states is 
given by 



k 



UJ 



-E,(k)+iT Lu + E,(k)+iT 



r is the quasiparticle scattering rate which origins from 
the lifetime effects, stoichiometry variations, noise smear- 



ing, etc 



16.17 



It can be also taken as a free parameter 



associated with the energy resolution in the tunneling 
experiment if the scattering rate is smaller than the ex- 
perimental resolution. 

Since the low-energy physics is governed by excitations 
near the Fermi energy, we assume that the kinetic en- 
ergy depends only on the absolute value of the momen- 
tum, i.e., e(k) = e{k). The anisotropic d-wave gap func- 
tion and c-axis coupling t±(k) can also be simplified as 
Ai(k) = Ai^o cos(2(?!)) and tj_{'k) = —tzCos'^{24>) in the 
vicinity of the Fermi surface. 

Near the Fermi energy, the normal density of states 
can be written as 



PAr(e) ~ PAr(O) + p'^r(0)e 



(9) 



up to the leading order approximation in e, where pjv(O) 
is the normal DOS at the Fermi energy. /ci^(O) is the 
linear coefficient of the DOS. p^(0) is finite if particle- 
hole symmetry is broken. A number of ARPES exper- 
iments have shown that there is a flat band at about 
200 meV below the Fermi energy in deeply underdoped 



cuprates 



18,19 



The presence of this flat band is an in- 



variation of the DOS around the Fermi surface can no 
longer be neglected as in conventional metals. Moreover, 
as revealed by the tunneling measurements, the tunneling 
conductance varies almost linearly with the applied bias 
around the Fermi energy in the normal stateii^ It sug- 
gests that particle-hole symmetry is broken in high-Tc 
cuprates. Thus, it is important to include the particle- 
hole asymmetric term in the analysis of tunneling mea- 
surement data. The linear approximation of DOS, de- 
fined by Eq. @, is valid if the energy of the Van Hove 
singularity is close to the Fermi energy but still much 
lower than the energy range we are interested in. 

With the above equations, it is straightforward to show 
that 

Pi{uj) ~ pAr(0)/i,i(w) + p^(0)sgn(tj)/^_2(w) 

±tMmM, (10) 



where Ii^i{uj) are 



27r2 



de- 



de- 



, (11) 



, (12) 



9 , (13) 

tj|-a) +r2 



dication of Van Hove singularity and suggests that the 



and rii = y^Afp cos2(20) + . The flrst term at 

the right-hand side of Eg. HlUI) has the largest contribu- 
tion to the DOS. The supperconducting coherence peaks 
are located at the gap edge, namely at tOi — iA^^o- 
In the limit F 0"*", liA^) oc |a;| near the Fermi 
energy, thus the low energy DOS of quasiparticles is 
linear pi(uj) ~ pAr(O) |w| /A^^o- These special d-wave 
characters have already been observed in the tunneling 
measurementsSiSiSiiSiiiiiSii^. 

Particle-hole symmetry is broken by the second term 
in (|10|l . The asymmetric DOS induced by this term has 
been observed in STM or other tunneling spectraiSiS 
Since liA^) is positive, a negative jO^(O) will enhance 
the DOS below the Fermi energy, but reduce that above 
the Fermi energy. 

The bilayer coupling appears in the third term of Eq. 
(|10|l . It leads to the difference in the DOS of the bonding 
and antibonding bands. This difference is proportional 
to both the c-axis hopping integral tz and p'jqiQ)- For a 
system with p^(0) < 0, the bonding band DOS pi(w) 
is reduced and the antibonding band DOS P2{'^) is en- 
hanced. 

Because of the c-axis coupling, the coherence peaks 
of the bonding and antibonding bands are separated. 
The maxima of the energy gap for both bonding and 
antibonding bands are located near the four points 
(±7r, 0), (0, ±7r). Therefore, the largest difference in the 
energy gap, 5^ = — ^2.,o\^ also occurs near these 

positions, where A^^o — ^(ki,F) is the gap maximum 
of the i'th band at the Fermi momentum k^ along the 
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Fig 1: Density of states in a bilayer compound. The param- 
eters used are Ai,o ~ 30 meV, A2,o = 33 meV, tz — 50 
meV, r = 0.25 meV, piv(O) = 1 eV-\ (a) p'jv(O) = 0, (b) 
p'jv(O) = -8 eV-^ 
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Fig 2: Same as for Fig. 1 but with Ai,o = 33 meV, A2,o = 30 
meV. 



anti- nodal direction. From the ARPES data in Ref.[4], 
we estimate the difference between the anti-nodal Fermi 
momentum to be \\i.i.F — k2,F| — 0.127r and the differ- 
ence of the largest energy gaps to be Sg ~ 0.1 Aq. In OD 
Bi2212 compounds, Aq ~ 30 meV, thus Sg ~ 3 meV. 

Fig. 1 shows the DOS for a system with Ai.o = 30 
meV, A2.0 = 33 meV, and PAr(O) — 1 eV^^. The results 
for both p^(0) ~ and /O^(0) = —8 eV^^ are shown. 
In the case Pjv(O) — 0; two coherence peaks appear due 
to the bilayer splitting. These two peaks, at 30 meV 
and 33 meV, come from the bonding and antibonding 
bands, respectively. When p'^(O) = —8 eV~^, the DOS 
becomes non-symmetric, as expected. This asymmetric 
feature exists in both the normal and superconducting 
states, consistent with the reported dataiiSiH However, 



the peak of the bonding band is strongly suppressed. 

For comparison. Fig. 2 shows the DOS for a system 
with Ai o — 33 meV and A2,o — 30 meV. The asymmet- 
ric behavior shown in Fig. 2 (b) is similar to the case for 
Fig. 1 (b), but the higher coherence peak is located at 
lower ene rgy. 
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Fig 3: Density of states with (a) F = 0.5 meV (b) F = 0.75 
meV. Other parameters used are the same as for Fig. 1 (b) 



Compared with the case Pjv(O) = 0, the particle- hole 
asymmetry suppresses strongly the lower coherence peak. 
If the quasiparticle scattering rate F is large or the energy 
resolution is not high enough, it is certainly difficult to 
resolve this double-coherence-peak structure in the tun- 
neling spectra. This is explicitely illustrated in Fig. 3. 
With increasing F, the double-peak structure disappears 
gradually and the lower coherence peak becomes indis- 
tinguishable from the background. Therefore, in order 
to observe this bilayer splitting in tunneling spectra, ex- 
perimental measurements with high quality single crystal 
and high energy resolution are desired. 

In conclusion, the superconducting coherence peaks 
are separated in bilayer high-Tc superconductors and can 
be used to probe the bilayer coupling effect with tunnel- 
ing measurements. The particle-hole asymmetry in the 
DOS enhances one of the coherence peaks, but reduces 
another one. If the life-time of quasiparticles is very short 
or the energy resolution in tunneling measurements is not 
high enough, the lower coherence peak is difficult to be 
resolved from the conductance background. This gives 
a qualitative account why the bilayer splitting has not 
been unambiguously observed in tunneling spectra. 
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